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Abstract 

Generalized Bose-Einstein (BE) and Fermi-Dirac (ED) distributions in 
nonextensive quantum statistics have been discussed by the maximum-entropy 
method (MEM) with the optimum Lagrange multipher based on the exact in- 
tegral representation [Rajagopal, Mendes, and Lenzi, Phys. Rev. Lett. 80, 
3907 (1998)]. It has been shown that the {q — 1) expansion in the exact ap- 
proach agrees with the result obtained by the asymptotic approach valid for 
0{q — l). Model calculations have been made with a uniform density of states 
for electrons and with the Debye model for phonons. Based on the result 
of the exact approach, we have proposed the interpolation approximation to 
the generalized distributions, which yields results in agreement with the exact 
approach within 0{q — 1) and in high- and low-temperature limits. By using 
the four methods of the exact, interpolation, factorization and superstatisti- 
cal approaches, we have calculated coefficients in the generalized Sommerfeld 
expansion, and electronic and phonon specific heats at low temperatures. A 
comparison among the four methods has shown that the interpolation approx- 
imation is potentially useful in the nonextensive quantum statistics. Supple- 
mentary discussions have been made on the {q—1) expansion of the generalized 
distributions based on the exact approach with the use of the un-normalized 
MEM, whose results also agree with those of the asymptotic approach. 
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1 INTRODUCTION 

In the last decade, many studies have been made for the nonextensive statistics 
[1] in which the generahzed entropy (the TsaUis entropy) is introduced (for a re- 
cent review, see [2]). The Tsalhs entropy is a one-parameter generahzation of the 
Boltzmann-Gibbs entropy with the entropic index q: the Tsalhs entropy in the limit 
of g = 1.0 reduces to the Boltzmann-Gibbs entropy. The optimum probability dis- 
tribution or density matrix is obtained with the maximum entropy method (MEM) 
for the Tsallis entropy with some constraints. At the moment, there are four pos- 
sible MEMs: original method [T], un- normalized method |3j, normalized method 
in, and the optimal Lagrange multiplier (OLM) method [5]. The four methods are 
equivalent in the sense that distributions derived in them are easily transformed 
each other [6]. A comparison among the four MEMs is made in Ref. [2]. The 
nonextensive statistics has been successfully applied to a wide class of subjects in 
physics, chemistry, information science, biology and economics [7]. 

One of alternative approaches to the nonextensive statistics besides the MEM 
is the superstatistics [H [9] (for a recent review, see [IQ]). In the superstatistics, it 
is assumed that locally the equilibrium state is described by the Boltzmann-Gibbs 
statistics and that their global properties may be expressed by a superposition over 
the intensive parameter (i.e., the inverse temperature) [H]-[in]- It is, however, not 
clear how to obtain the mixing probability distribution of fluctuating parameter from 
first principles. This problem is currently controversial and some attempts to this 
direction have been proposed [II]-[15]. The concept of the superstatistics has been 
applied to many kinds of subjects such as hydrodynamic turbulence |l6l [TTl [18], 
cosmic ray [I^ and solar flares [20]. 

The nonextensive statistics has been applied to both classical and quantum sys- 
tems. In this paper, we pay attention to quantum nonextensive systems. The 
generalized Bose-Einstein and Fermi-Dirac distributions in nonextensive systems 
(referred to as g-BED and g-FDD hereafter) have been discussed by the three meth- 
ods, (i) The asymptotic approximation (AA) was proposed by Tsallis, Sa Barreto 
and Loh [21] who derived the expression for the canonical partition function valid 
for \q — l\/kBT ^0. It has been applied to the black-body radiation [21], early 
universe |2H 122] and the Bose-Einstein condensation [21] [23] . (ii) The factorization 
approximation (FA) was proposed by Biiyiikkilic, Demirhan and Giilec [21] to eval- 
uate the grand canonical partition function. The FA was criticized in f2^ [SB] , but 
supported in [27], related discussion being given in Sec. 4. The simple expressions 
for g-BED and g-FDD in the FA have been adopted in many applications such as 
the black-body radiation [23l [28], [29l [30], early universe [STJ [32], the Bose-Einstein 
condensation [33]- [39], metals [iQ], superconductivity [HI [12], spin systems [l3]-[i8] 
and metallic ferromagnets [IH]. (iii) The exact approach (EA) was developed by 
Rajagopal, Mendes and Lenzi [501121] who derived the formally exact integral repre- 
sentation for the grand canonical partition function of nonextensive systems which is 
expressed in terms of the Boltzmann-Gibbs counterpart. The integral representation 
approach originated from the Hilhorst formula [52J. Because an actual evaluation 
of a given integral is generally difficult, it may be performed in an approximate way 
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[sot [51] or in the limited cases [53] . The vahdity of the EA is discussed in [511 EH] • 
The EA has been apphed to nonextensive quantum systems such as black body 
radiation [56l [57] and the Bose-Einstein condensation [501 EI] • 

We believe that it is important and valuable to pursue the EA despite its diffi- 
culty. It is the purpose of the present study to apply the EA [HUl [SI] to calculations 
of the generalized distributions of g-BED and g-FDD. The grand canonical partition 
function of the nonextensive systems is derived with the use of the OLM scheme 
in the MEM [5]. Self-consistent equations for averages of the number of particles 
and energy and the grand-canonical partition function are exactly expressed by the 
integral representation [501 [51] • The integral representation for q > 1.0 in the EA 
is expressed as an integral along the real axis, while that for g < 1.0 is expressed 
as the contour integral in the complex plane [50l [STl [53] . We have shown that the 
{q — 1) expansion by the EA agrees with the result derived by the AA. For q > 1.0, 
the self-consistent equations have been numerically solved with the band model for 
electrons and the Debye model for phonon. 

It is rather difficult and tedious to obtain the generalized distributions in the EA 
because they need the self-consistent calculation of averages of number of particles 
and energy. Based on the exact result obtained, we have proposed the interpolation 
approximation (lA) to g-BED and g-FDD, which do not need the self-consistently 
determined quantities and whose results are in agreement with those of the EA 
within 0{q — l) and in high- and low-temperature limits. We may obtain the simple 
analytic expressions of the g-BED and g-FDD. 

The paper is organized as follows. In Sec. 2, the exact integral representation is 
derived with the OLM-MEM after Ref. [501 [SH [53] . We have discussed the (g - 1) 
expansion of physical quantities, using the EA and AA. Numerical calculations are 
performed for electron and phonon models, for which we present the g-BED and 
g-FDD with the temperature-dependent energy. In Sec. 3, we propose the lA, 
by which analytical expressions for g-BED and g-FDD are obtained. In Sec. 4, 
a comparison is made between the generalized distributions calculated by the four 
methods of the EA, lA, FA [21] and the superstatistical approximation (SA). A 
controversy on the validity of the FA [23] is discussed. With the use of the four 
methods, the generalized Sommerfeld expansion, and low-temperature electronic 
and phonon specific heats are calculated. Sec. 5 is devoted to our conclusion. In 
Appendix A, we present a study of the EA and AA with the un-normalized MEM 
[3l[2T], calculating the (g — 1) expansion of the g-BED and g-FDD. Supplementary 
discussions on the lA are presented in Appendix B. 

2 EXACT APPROACH 
2.1 MEM by OLM 

We will study nonextensive quantum systems described by the hamiltonian H. We 
have obtained the optimum density matrix of p, applying the OLM-MEM to the 
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Tsallis entropy given by [5i i6j 
with the constraints: 

Trpg = 1, 
Tr{plN} = c,N„ 
Tr{plH] = c,E,, 
Cq = Trpl, 

where Tr stands for the trace, ks is the Bohzmann constant, and Eg and Ng denote 
the expectation values of the hamiltonian H and the number operator N, respec- 
tively. The OLM-MEM yields [3 E] 

Pg = -^[l + iq-l)P{H-pN-Eg + pNg)]Th;^ (1) 

Xg = Tr{[l + iq-l)/3iH-fiN-Eg + pNg)]Tk}, (2) 
Ng = -LTr{[l + iq-l)P{H-pN-Eg + pNg)]T^. N}, (3) 

Eg = -LTr{[l + {q-l)P{H-pN-Eg + pNg)]^^ H}, (4) 

where /? and p denote the Lagrange multipliers. In deriving Eqs. ([I])-(|4]), we have 
employed the relation: 

^g 

Lagrange multipliers of (3 and p are identified as the inverse physical temperature 
(/3 = l/Zc^T) and the chemical potential (fermi level), respectively. OE]. 

2.2 Exact integral representation 
2.2.1 Caseofg>l 

In the case of g > LO, we adopt the formula for the gamma function r(s): 

= -— / n^-^e"^" du for s > 0. (5) 

r(s) Jo 

With s = l/(g - 1) [or s = q/{q - 1)] and x = 1 + (g - l)f3{H - pN) in Eq. ©, we 
may express Eqs. ©-dl]) by [50} [51] 



Xq Jo V ? - 1 



(6) 
(7) 



4 



with 
where 

k 

Qiiu) = i-V ln[l Te-"(^^-^)], (10) 

k 

Ni{u) = 5^/i(e,,M), (11) 

k 

E,{u) = J2^kfiiek,u), (12) 

k 

/^(^'-) = ^MT^V ^''^ 

G{u;a,b) = -^u^-'e-'\ (14) 
r(a) 

The upper (lower) sign in Eqs. ([9]), (fTO!) and (fT3!) denotes boson (fermion) case, 
and Si(m), f2i(M), Ni{u), Ei{u) and fi{e,u) express the physical quantities for q = 
1.0. Equations ([S])-(IH]) show that physical quantities in nonextensive systems are 
expressed as a superposition of those for q = 1.0. 

Although Eqs. ([6])- ([8]) are formally exact expressions, they have a problem when 
we perform numerical calculations. The gamma distribution of G[u] l/(g — 1) + £, 1] 
{i = 0, 1) in Eqs. ([n])-(IHD has the maximum at Umax, and average and variance given 
by 

Umax — 7 TT ~\~ ^ 1; (1^) 

w- 1) 

= (16) 

{u\-{u)l = (^ + ^- (17) 

Equation f|T5l) shows that the gamma distribution in Eqs. ([6])- ([8]) has the maximum 
at Umax = l/{q — 1) ^ oo while the contribution from Si[(g — l)/3t] is dominant at 
t ~ because its argument becomes (g — l)f3t 0. Then numerical calculations 
using Eqs. (I6])-(IHD are very difficult. 

In order to overcome this difficulty, we have adopted a change of variable: (g — 
l)j3u — > M in Eq. ([6])-([8]) to obtain alternative expressions given by 
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with 

The gamma distribution of G{u; + £, ^/?) for £ = 0, 1 in Eqs. ([I8D-([20D has 
the maximum at Umax, and average, mean square and variance given by 

Umax = + (21) 

{u)u = + (22) 

{u')^ = [l + (q-l)i][l + (q-l)(i+l)]p^, (23) 

{n')u-{u)l = iq-l)[l + {q-m(3'. (24) 

Equation (12T!) shows that the gamma distribution has the maximum at Umax = (3 
in the hmit of q —* 1.0, and an integration over u in Eqs. ( IT8|) - (!20|) may be easily 
performed. Indeed, in the case of g > 1.0 discussed above, the gamma distribution 
in Eqs. (HHD-d^DD becomes 

V'g-1 '{q-l)PJ V27r(g-l)/32 

5(u - (3) for (g - ^ 0. (26) 

Although expressions given by Eqs. ([6])- ([8]) are mathematically equivalent to those 
given by Eqs. ffT51) - fl2Ul) . the latter expressions are more suitable than the former 
ones for numerical calculations. 

2.2.2 Case of g < 1 

In the case of g < 1.0, we adopt the formula given by 

x' = — r(s + l) /(-t)-^-^e-^*dt forSRs>0, (27) 
2vr Jc 

where a contour integral is performed over the Hankel path C in the complex plane. 
With s = 1/(1 - g) [or s = g/(l - g)] and x = 1 + (g - 1)/3{H - fiN) in Eq. 
we obtain |50l 



= ^^^(t;Y3^,l)e-(i-^)^*(^'-^^')Si[-(l-g)/5t]iVi[-(l-g^^^^ 

(28) 

= ^^^(t;Y^,l)e-(^-^)^*(^^-'^^^)Hi[-(l-g)/3t]Ei[-(l-g) 

(29) 



with 



I f ..f 1 

H{t;a,b) = r{a + l)b-^ {-t)-''-^e-'", (31) 



X, = -y^if(^t;^,lJe-(i-^)^*(^'-'^^^)Si[-(l-g)/3t]dt, (30) 



6 



where Ni{u) Ei{u) and fi{e,u) are given by Eqs. (!9|)- (fT2|) with complex u. 

In the case of g < 1.0, Ng, Eg and Xg given by Eqs. fl28|) - (!30l) are expressed 
by an integral along the Hankel contour path C in the complex plane. The Hankel 
path may be modified to the Bromwich contour which is parallel to the imaginary 
axis from c — iootoc + ioo (c> 0) [Sni The Bromwich contour is usually 
understood as counting the contributions from the residues of all poles located in 
the left-side of 3? z < c of the complex plane z, when the integrand is expressed 
by simple analytic functions. If the integrand is not expressed by simple analytic 
functions, we have to evaluate it by numerical methods. Unfortunately, we have not 
succeeded in evaluating Eqs. (I28l)-( l30|) with the sufficient accuracy. It is not easy 
to numerically evaluate the integral along the Hankel or Bromwich contour, which 
is required to be appropriately deformed for actual numerical calculations [581 ES] • 
This subject has a long history and it is still active in the field of the numerical 
methods for the inverse Laplace transformation [58] and for the Gamma functions 



It is worthwhile to remark that for a bose gas model with the density of states 
of p(e) = Ae'', we obtain (with /x = 0) [521 [56l [57] 

"Ar(r + l)C(r + 2) 



^i[u) = exp 

Ar(r + l)C(r + 1) 



EAu) 



Ar(r + 2)C(r + 2) 



where r = 1/2 for an ideal bose gas, r = 2 for a harmonic oscillator, A denotes a 
relevant factor and C,{z) stands for the Riemann zeta function. With a repeated use of 
Eq. (P7|) . Nq, Eg and Xg may be expressed as sums of gamma functions [S21 ESI EZ]- 
Unfortunately, such a sophisticated method cannot be necessarily applied to any 
models like a fermi gas. 

With a change of variable of (1 - q)l3{-t) {-t) in Eqs. (128D-([30D after the 
case of g > 1, they are given by 

(32) 



(33) 



with 



Average and mean square over H(t, — £, (jr^) for i = 0,1 are given by 

{i-t))t = [l-(l-gm (35) 
{{-t)% = [l-{l-q)i][q-{l~qW. (36) 
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Equations (!32l) - (!MI) are useful in making the — 1) expansion, as will be discussed 
in the following. 



2.3 The [q — 1) expansion 
2.3.1 The exact approach 

We will consider the (g — 1) expansion of the expectation value of an operator O in 
the EA. By using Eqs. ([18]) and ([32]), we obtain 



X 



Tr {[1 - {1 ~ q)(3k]^^ 0}, 



— I Giu; 



Jo 

i 



27rX 



H t 



q JC 



g-l'(g-l)/3 
q 1 



l-g'(l-g)/3 



(37) 

Yi{u)Oi{u)du forg>l, (38) 
\Yi{-t)Oi{-t)dt forg<l, (39) 



with 



Oi(m) 
Y,{u) 

k 



Tr{e-''^ 0} 
Yiiu) ' 

H -jjN -Eq + jjNq, 



(40) 

(41) 
(42) 



where Xg is given by Eq. f[20[) for g > 1 and by Eq. f[34[) for g < 1. It is noted that 
Yi{u) includes the self-consistently calculated Ng and Eg. 

We first consider the case of g > 1 for which the integral including an arbitrary 
function W{u) is assumed to be given by 



J 



oo / ^ 

G I u: + 







g-1 '(g-l)/3. 



W{u) du 



for ^ = 0, 1. 



(43) 



Since G{u; +i, (~ij^) has the maximum around u = P as mentioned before [Eq. 
( [2T[) ]. W{u) may be expanded as 



dW 1 ^d'^W 
W{u) = W{f3) + {u-f3)^ + -{u-py^ + --. (44) 

Substituting Eq. ([H]) to Eq. ([151) and using the relations given by Eqs. ( [221) and 
), we obtain J in a series of (g — 1) as 



dW 1 o d'^W 



W{P) + (g - 1) 



+ 



(45) 

for g ~ 1.0. (46) 
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Next we consider the case of g < 1 for which a similar integral along the Hankel 
path C is given by 



c 



l-q 



W{-t) dt for £ = 0, 1. (47) 



By expanding W{—t) aX —t = (3 and using the relations for averages given by Eqs. 
(1351) and (136|) . we obtain the same expression for J as Eq. (H6l) . which is then valid 
both for q < 1.0 and q > 1.0. 

For W{u) = Yi{u) and W{u) = Yi{u)Oi{u) in Eq. (iEl), we obtain 



df3^ 



+ 



(48) 
(49) 



Note that the 0{{q — l)j3) term in Eq. (HHl) vanishes because £ = in Eq. (j46] 
Substituting the relations given by 



dYi 

dOi 
df3^ 



-{K),Y,, 

(ir)i(6)i-(K6)i, 



to Eqs. 



and din]), we finally obtain the 0{q — 1) expansion of Og given by 



Wq — 



Oi + (1 - q) (/3(i^0)i + l/3'[{K'h{dh - {K'Oh] I + 



(50) 
(51) 
(52) 
(53) 

(54) 



2.3.2 The asymptotic approach 

On the other hand, we may adopt the AA [21] to obtain Og given by Eq. fl371) valid 
for 0{q — 1). By using the relation: ~ e^[l — (1 — q)x^/2 + ■■] in Eqs. ([2]) and 
( 1571) . we may expand and Og up to 0(g — 1) as 



On 



l--{l~q)f3'{K\ + 



^Tr{[l - (1 - q)pk]-' [1 - (1 - q)PK]^ O}, 



X, 



Tr{e-'^^[l + {l-q)(3K] 



l-^{l-q)P^K^ 



0} + -, 



Oi + (1 - g) ( /?(A-6)i + -P' [{K')i{d), - {k^O), 



+ 



(55) 
(56) 
(57) 
(58) 



Equation (!58|l agrees with Eq. (15^ obtained by the EA within 0{q—l). In Appendix 
A, we have shown that the same equivalence holds between the AA and EA with 
the un-normahzed MEM p[2T]. 



2.4 Generalized distributions 
2.4.1 The 0{q — 1) expansion 

Equations for Ng and Eg given by Eqs. ( fTSi) . ( fT9l) . ( !32|) and ( !33l) may be expressed 



as 



= $^/,(efc,/5)= / fg{e,P)pie)de, 
k 

Eq = ^fqiek,P)ek= fg{e,P)ep{e)de, 

7, *^ 



(59) 
(60) 



where fg{e,P) [= /g(e)] signifies the generalized distributions, g-BED and g-FDD, 
given by 

for g > 1, (61) 



H t 



1 



2'KXgJc yi-q' {l-q)(3 
with the density of states p(e) given by 



p(<) 



ri(-t)/i(e, -t) 

for g < 1, (62) 

(63) 



In order to examine the (g — 1) expansion of the generalized distributions, we set 
O = nfc in Eq. (15^ where fik denotes the number operator of the state k. A simple 
calculation leads to the 0{q — 1) expansion of the generalized distribution given by 



fg{e,P) = f,{e,P) + {q-l) 
= /i(e,/3) + (g-l) 



^ d(3 r d(3^ 



+ ■ 



(64) 
(65) 



In deriving Eq. ( l65l) . we have employed the relation: ((9Yi/9/5)/Y'i(/3) = —{H — 
/iA^)i + {Ng — iiNg) ~ 0(g — 1). In Appendix A, we have made a similar analysis 
with the un-normalized MEM, showing that Eq. fl65l) is consistent with Eq. (]A39I) 
which agrees with the result in the AA [21] . 
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2.4.2 Properties of the generalized distribution 



We will examine some limiting cases of the generalized distribution given by Eqs. 
(|6T|) and ([62]). 

(1) In the limit of g ^ 1.0, Eq. leads to 



f,{e,P) = /i(e,/3). (66) 
(2) In the zero-temperature limit of /? — oo, the g-FDD becomes 



/,(e,r = 0) = e(/i-e) = /i(e,T = 0), (67) 

where 0(x) stands for the Heaviside function. Equation (1671) implies that the 
ground-state FD distribution is not modified by the nonextensivity. 
(3) In the high-temperature limit of /3 — *• 0.0, where ili ~ — (1//5) '^f, e~^^^''~'^^ with 

ln(l ± x) ~ =i=x for small x, we obtain (/x = 0.0) 

f,{e,(3^0) cc [l + iq-l)f3{e-E,)]^.-' = [e~^'~^-^'>]\ (68) 

expressing the g-exponential function defined by 

Cg = expg(a;) = [1 + (1 - g)x]i^ for 1 + (1 - g)a; > 0, (69) 
= for 1 + (1 - q)x < 0, (70) 

with the cut-off properties. Equation (!68l) corresponds to the escort distribution, 

W = ^^[ef^'-'-^r, (71) 

Cq 

with the g-exponential distribution Pq{e) given by 

Pq{e) = e~^^^~^\ (72) 

Equations (I6T1) and (!62l) shows that the e dependence of fq{e, (3) arises from that 
of /i(e, /?). In particular, the g-FDD preserves the same e symmetry as /i(e, /5): 

(a) /,(e,/?) = l/2fore = /i, 

(b) fq{e,P) has the anti-symmetry: 

fq{-Se + ^, /5) - 1 = 1 - fg{Se + fi, P) for Se > 0, 

(c) dfq{e,/3)/de is symmetric with respect to e = /i. 
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2.5 Numerical calculations 



2.5.1 Model for electrons 

For model calculations of electron systems, we employ a uniform density of state 
given by 

p(e) = (l/2W^)e(W^-|e|), (73) 

where W denotes a half of the total band width. We have performed numerical 
calculations of Eg and fi for q > 1.0 as a function of T for a given number of 
particles of N and the density of states p(e) . We may obtain analytical expressions 
for Ni{u) and Ei{u) which are necessary for our numerical calculations. By 

using Eq. ([73]) for Eqs. ©-(US), we obtain (with W = 1.0) 

Si(u) = e-"^i("), 

Qi{u) = + e-"^^-'^)] - ln[l + e-"(^+^)] + ln[l + e"(^+^)] - ln[l + e"^^-'^)]} 



1 

2^ 



{L22(-e-"(^+^))-L22(-e"(^-^))}, 
1 



N^{u) = l + ^[ln(l + e-"(^+'^))-ln(l + e"(^-'^))], 
E^{u) = -i-[ln(l + e-"(i+^))+ln(l + e"(i-'^))] 

+ -i^[Lz2(-e-"(i+'^)) -L22(-e«(i-'^))], 
where Lin{z) denotes the nth polylogarithmic function defined by 

oo p, 
k=l 

We adopt = 0.5, for which fi = 0.0 independent of the temperature because 
of the adopted uniform density of states given by Eq. (1751) . The temperature 
dependence of Eg calculated self-consistently from Eqs. f|T8l) - fl20l) . is shown in Fig. [1] 
whose inset shows the enlarged plot for low temperatures {kBT/W < 0.1). We note 
that Eg at low temperatures is larger for larger q although this trend is reversed at 
higher temperatures (fc^T > 0.3). 

The calculated g-FDDs fg{e) for various q values for ksT/W = 0.1 are shown 
in Figs. [2] (a) and [2] (b) whose ordinates are in the linear and logarithmic scales, 
respectively. It is shown that with more increasing q from unity, /g(e) at e ^ /i 
has a longer tail. The properties of fg{e) are more clearly seen in its derivative of 
—dfg{e)/de, which is plotted in Fig. [3] with the logarithmic ordinate. We note that 
—dfg{e)/de is symmetric with respect of e = /i. With increasing q above unity, 
—dfg{e)/de has a longer tail. Dotted and sohd curves for q < 1.0 in Figs. [2] and [3] 
will be discussed in Sec. 3.3. 
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2.5.2 The Debye model for phonons 

We adopt the Debye model whose phonon density of states is given by 

p{uj) = AuP' for < cij < ujd, 



(74) 



where A = 9Na/w^, Na denotes the number of atoms, u the phonon frequency and 
ujd the Debye cutoff frequency. By using Eq. (!7i|) to Eqs. (!9!)- (fT2l) . we may obtain 
(with uj£, = 1.0 and /i = 0), 



'U Hi (ji) 



A 



180m 



47r^ 



+ 15u — 60 ln(l — e") + 60 ln(l — cosh u + sinh u) 



A 
u 



4 [u^ L22(e") - 2u Lz3(e") + 2L^4(e")], 



A 



-^[m^ - 3M^ln(l - e") - 6MLi2(e") + 6Lz3(e") - 6 C(3)] 
In(l-e^) 3Li2(e") 6Li3(e") 6Li4(e") 1 



TT 
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We have performed numerical calculations with the Debye model for q > 1.0. 
The temperature dependence of self-consistently calculated Eg is shown in Fig. H] 
where inset shows the enlarged plots for low temperatures {T/Tj:, < 0.5). We note 
that Eg at low temperatures is larger for larger q. 

The calculated g-BEDs fg{e) for various q values for T/Td = 0.01 are shown in 
Fig. whose ordinate is in the logarithmic scale: they are indistinguishable in the 
linear scale. It is shown that with more increasing q, fg{e) at e ^ /i has a longer 
tail. Dotted and solid curves for g < 1.0 will be discussed in Sec. 3.3. 



3 THE INTERPOLATION APPROXIMATION 

3.1 Analytic expressions of the generalized distributions 

In the preceding Sec. 2, we have discussed the generalized distributions based on 
the exact representation given by Eqs. (16T|) and (1621) . It is, however, difficult to 
calculate them because they need self-consistent calculations of Ng and Eg. If we 
assume 

e"(^^-'^^^)Si(M) = 1, (75) 
in Eqs. ( !6T!) and (|62|) . we obtain the approximate generalized distributions given by 
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where G{u; a, h) and H{t] a, h) are given by Eqs. and (!3T|) . respectively. Equa- 
tions (!76l) and (!77|) are referred to as the interpolation approximation (lA) in this 
paper because they have the important interpolating character, as will be shown 
shortly (Sec. 3.2). Note that calculations of /^^(e, /?) by Eqs. (175]) and (I77|) do not 
require Nq and Eg. Equation (1761) may be regarded as a kind of the SA. 

One of advantages of the lA is that we can obtain the simple analytic expressions 
for the g-BED and g-FDD as follows. 
(1) g-BED 

We first expand the Bose- Einstein distribution /i(e, /3) as 



oo 



/i(e,/3) = ^ e-("+^)" forx>0, (7J 



n=0 



where x = P{e — fi). Substituting Eq. flTHl) to Eqs. fl76l) and (177|) . and employing 
Eq. ([5]) and fl27|) . we obtain the g-BED in the lA given by 



C{e,f3) = $^[e-("+^)T for0<g<3, (79) 

n=0 

_ 9 

7 1 



1 


<3 


-c( 


(g — l)x 



^g — 1 ' (g — l)x 
where ({z, a) denotes the Hurwitz zeta function: 



1 for 1 < g < 3, (80) 



((z, a) = y — — = — — / dt for 3? z > 1. 



,=0 (k + ay T{z)Jo 1-e- 



It derivative is given by 



dx 



g(n + l)[e-("+^)^](2'?-i) for < g < 3. (81) 



n=0 



We may easily realize that fq{e,P) in Eq. fl79l) reduces to fi{e, (3) in the limit of 
g — > 1.0 where — > e^. 

(2) g-FDD 

The Fermi-Dirac distribution /i(e, /?) may be expanded as 



h{e,(5) = ^(-l)«e-("+i)- forx>0, (82) 

for x = 0, (83) 
^ (-l)"e~''l^'l fora;<0, (84) 



n.=0 
1 

2 



n=0 
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where x = /5(e - /x). Substituting Eqs. ^-([HID to Eqs. ([76]) and ([77D, and 
employing Eq. ([5]) and (!27|) . we obtain the g-FDD in the lA given by 



with 



F,(x) 



n=0 



1 
2 

l-^.(kl) 



-1)" [e-("+i)^l'? 



_2{q-l)x 



q-l 



{C 



for X > 0, 

for a; = 0, 
for a; < 0, 

for < g < 3, 
1 



^5) 

^6) 
^7) 



-c 



q-l' 2{q-l)x 



1 

g-l'2(g-l)x ^ 2 



+ 1 } for 1 < g < 3. 



It derivative is given by 



dx 



q{n + l)[e- 



(n+l) |xh(2<7-l) 



for < g < 3, 



^9) 



(90) 



n=0 



which is symmetric with respect to x = 0. The g-FDD given by Eqs. (!85i) - (!88l) 
reduces to /i(e, /5) in the hmit of g — 1.0. 

We may obtain a useful expression of the g-FDD for |x| < 1 given by (see 
Appendix B.l) 



lA 



dx 



1 q , g(2g-l)(3g-2) 3 
X H X + 

2 4 48 

g , g(2g-l)(3g-2) , 
4+ 16 + 



for < g < 3. 



(91) 
(92) 



In the case of g < 1.0, summations over n in the g-BED and g-FDD [Eqs. (1791) 
and (1551) ] are terminated when the condition: n + 1 > 1/(1 — q)x is satisfied because 
of the cut-off properties of the g-exponential function given by Eq. fl70|) . Then the 
g-FDD for g < 1.0 has the cut-off properties given by 



0.0 for e - /i > 1/(1 - g)/5, 
1.0 for e-/i < -1/(1 -g)/3, 



(93) 
(94) 



while the g-BED has the cut-off properties given by Eq. fl93|) . These are the same 
as the g-exponential distribution Pq{e) given by Eq. fl72|) . 



3.2 Comparison with the exact approach 

From Eqs. (08]) and (gg]) with Yi{u) = 1.0, the g-BED and g-FDD for g c:^ 1.0 in the 
lA become 

r pti- 1 - 

cIAi 



C{e,f3) = Me, (3) + (q-l) 



9/1 1 ,d'f: 



de 



+ 



(95) 
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which is in agreement with those in the EA given by Eq. (!65|) within 0{q — 1). In 
the zero-temperature hmit, the g-FDD reduces to 




e,T = 0) = e(/i-e). 



(96) 



In the opposite high-temperature limit, the g-BED and g-FDD become 



//^(6,/5^0)(x[e-T. 



(97) 



Equations (!96|) and (!97|) agree with Eqs. (1671) and (!68|) . respectively, for the EA. Thus 
the generalized distributions in the lA have the interpolation properties, yielding 
results in agreement with those in the EA within 0{q — 1) and in high- and low- 
temperature limits. 

3.3 Numerical calculations 

Numerical calculations of f^^{e,(3) [= fq^{e)\ have been performed. Results of the 
FDD of /f^(e) in the EA for q > 1.0 and ksT/W = 1.0 are shown in Fig. El 
With more increasing q, the distributions have longer tails, as shown in Fig. [2] for 
ksT/W = 0.1. The result in the lA is in good agreement with the EA because 
the ratio defined by A = /g^(e)//g^^(e) is 0.97 < A < 1.01 for -10 < e < 10 as 
shown in the inset. The e dependence of the BED of /^^(e) in the EA for q > 1.0 
and T/Tu = 0.1 is plotted in Fig. [7] which shows similar behavior to those for 
T/Td = 0.01 shown in Fig. [6l Its inset shows that the ratio of A is 0.7 < A < 1.0 
for 1.0 < q < 1.2. These calculations justify, to some extent, the distribution in the 
lA given by Eqs. ([80]), ([85])-([87D and (1891) . 

We have calculated the g-BED and g-FDD also for q < 1.0, by using Eqs. ( l79l) . 
dHSD-dHHD. Dotted and solid curves in Fig. El show the g-FDD of //^(e) for q = 0.9 
and q = 0.8, respectively. Their derivatives of —df^"^{e)/de for q = 0.9 and q = 0.8 
are plotted by the dotted and solid curves, respectively, in Fig. [31 Dotted and solid 
curves in Fig. [51 show the g-BED of /^'^(e) for q = 0.9 and q = 0.8, respectively. 
With more decreasing q from unity, the curvature of fq{e) in both g-BED and g-FDD 
become more significant. The cut-off properties in the g-FDD and g-BED given by 
Eqs. (1^ and (IMj) are realized in Figs. [21 and [3 We expect that fq^{e) in the 
case of g < 1.0 is a good approximation of the g-BED and g-FDD as in the case of 



4.1 Comparison with previous studies 

It is interesting to compare our results to those previously obtained with some 
approximations. 

(A) The factorization approximation 



q > 1.0. 



4 DISCUSSION 
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Biiyiikkilic, Demirhan and Giilec [2^ derived the g-BED and g-FDD given by 



adopting the FA given by 

AT 



Q = [l-(l-g)^x„] — , (99) 

n=l 

N 



n=l 



to evaluate the grand canonical partition function, the upper (lower) sign in Eq. 
f l98p being applied to boson (fermion). 

It is noted that if we assume the factorization approximation: [cq — 
{e-y[ie-'')X in /^^(e) [Eqs. (CHI) and ([HHD], we obtain 

which is similar to Eq. (!98|l [ill [55] . 

(B) T/ie super statistical approximation 

In the SA, the generalized distribution is expressed as a superposition of /i(e) 
[HI [9], 



/f^(e,/3) = / ^) /i(e,w)ciw, (102) 



which is similar to but different from (e, /3) given by Eq. (1751) . Recently the 
g-FDD equivalent to Eq. fl^JS]) is obtained by employing the SA in a different way 
[19]. 

The properties of the generalized distributions of the EA, lA, FA and SA in the 
limits ofg— i>1.0, /5— i>cxD and /? — 0.0 are compared in Table 1. The result of 
the lA agrees with that of the EA within 0{q — 1) as mentioned before. However, 
the 0{q — 1) contributions in the FA and SA are different from that in the EA. 
In the zero-temperature limit, all the g-FDDs reduce to B(yU — e). In the opposite 
high-temperature limit, the generalized distributions in the FA and SA reduce to 
e"^"^, while those in the EA and lA become [e"^'^]'' where the power index q arises 
from the escort probability in the OLM-MEM given by Eq. flTTl) [5\ 

Figure [HI shows g-BED for q = 1.1 and q = 1.2 calculated by the FA, SA and 
EA with the logarithmic ordinate. For a comparison, we show fq{e) for q = 1.0 by 
dashed curves. The difference among fq{e)^s of the three methods is clearly realized: 
tails in the g-BED of the FA and SA are overestimated. 

Figure [9] shows g-FDD for q = 1.1 and q = 1.2 calculated by the EA, FA and 
SA with the logarithmic ordinate (for more detailed fg^i^), see Fig. 1 of Ref. 
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Tails in the FA and SA have larger than that in the EA, as in the case of the g-BED 
shown in Fig. [HI 

Figures [TOT a) and 10(b) show the g-FDD and its derivative, respectively, calcu- 
lated in the lA and FA. For q = 0.9, fq^{e) at e < is much reduced than fg^{e). 
For g = 1.1, on the contrary, fg^{e) at e > yU is much increased than /^^(e). These 
lead to an overestimate of electron excitations across the fermi level /i in the FA. 
Furthermore —dfg^{e)/de in the FA is not symmetric with respect to e = in 
contrast to that in the lA. 

The FA was criticized in Refs. j25]j26] but justified in Ref. [27]. The dismissive 
study [25] was based on a simulation with N = 2. In contrast, the affirmative study 
[27] performed simulations with N = 10^ and 10^^. Lenzi, Mendes, da Silva and 
Malacarne p6| criticized the FA, applying the EA [50l [51] to independent harmonic 
oscillators with < 100. Our results are consistent with Refs. [251 126]. The FA 
given by Eq. fllOOl) has been explicitly or implicitly employed in many studies not 
only for quantum but also classical nonextensive systems. It would be necessary to 
examine the validity of these studies using the FA from the viewpoint of the exact 
representation ^50l [STj [60]. 

By using Eqs. ([5]) and ([27j) . we may rewrite Q in Eq. f[99l) as 

Q = [l-(l-g)xi]^®,--®Jl-(l-g)x^]^, (103) 

PCO / 1 1 \ ^ 

= / Giu; , TTe""'^"c^M for g > 1.0, (104) 

Jo V 9-V„=i 

= — f H (t;^—]Y\e'^^~'^^'''" dt for g< 1.0, (105) 

where ®g denotes the g-product defined by [62] 

x^gV = [x^-" + y^-'^ - 1]^^ . (106) 

Equations fll04p and (11051) are the integral representations of the g-product given by 
Eq. ( llOSp . The result of the FA in ( llOOp is derived if we may exchange the order of 
integral and product in Eqs. (11041) and (11051) . which is of course forbidden. 

4.2 The generalized Sommerfeld expansion 

We will investigate the generalized Sommerfeld expansion for an arbitrary function 
0(e) with the g-FDD of /^(e) given by 



(e)/,(e) de, (107) 

oo 

/)(e) de + ^c„,, {kBTr 0^"-'^/^), (108) 



n=l 



with 



cn, = - ^ /(^ - ^)"^ de. (109) 
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Substituting /q(e) in the EA given by Eq. (!65|) to Eq. (11091) . and using integrals by 
part, we obtain Cn,q for even n, 

= i + _!) + .. for even n, (110) 

= l + (g-l) + -- forn = 2, (111) 
= l + 6(g-l) + -- forn = 4, (112) 

while Cn^q = for odd n, where c„_i denotes the relevant expansion coefficient for 
q = 1.0:'c2,i = 7rV6 (=1.645) and C4,i = 7^^360 (=1.894) et. al. Equation ffTTOD 
shows that Cn,q is increased with increasing q. 

By using fg^{e) in the lA, we may obtain c„,g given by (for details, see Appendix 
B.2) 



c^^ r(^ + l-r2) 

n,q 1 ' 



C„,i (g-l)nr(-i^ + l) 

r(T^ + i) 



;i-g)«r(^ + i + n) 
1 



2-q 



(2-g)(3-2g)(4-3g) 



for even n, q > 1 , (113) 

for even n, q < 1, (114) 

for n = 2, (115) 

forn = 4. (116) 



It is easy to see that Eqs. (11151) and (I116p are in agreement with Eq. (Illip and 
(11121) . respectively, of the EA within 0(g — 1). 
A simple calculation using /^"^(e) leads to 

c^^ r(^ - n) 

Jhl = ^g-i for even n (q > 1), (117) 

(g-l)«r(^) ^' ^ ' 

for n = 2, 



(2-g)(3-2g) 
1 



(2-g)(3-2g)(4-3g)(5-4g) 



for = 4, 



which are similar to those given by Eqs. (11151) and (11161) . 

The Sommerfeld expansion coefficients in the FA may be calculated with the use 
of fq^{e) |l9]. A comparison among the 0{q — 1) contributions to Cn,q (n = 1 — 4) 
in the four methods of EA, lA, FA and SA is made in Table 2. The results of the 
lA coincide with those of the EA. The 0{q — 1) contributions to C2,q and C4 g in the 
SA are three and 5/3 times larger, respectively, than those in the EA. The 0{q — 1) 
contributions to C2,q and 04,^ in the FA are vanishing. It is noted that cf^^ 7^ and 
^3,q 7^ in contrast with the results of ci^q = c^^q = in the EA, lA and SA. This 
is due to a lack of the symmetry in —df^^{e)/de with respect to e = as shown in 

Fig. mh)- 



19 



Figure [TlTa) shows the q dependence of coefficients of Cn,q/ Cn,i for n = 2 and 4 
calculated by the four methods. Circles and squares express cf^ for n = 2 and 4, 
respectively, calculated by the EA for ksT/W = 0.1 (Fig. [1]). Solid curves express 
c^'^ in the lA. The coefficient for n = 2 (n = 4) in the lA is in good agreement with 
the result in the EA for 1.0 < q ^ 1.5 (1.0 < (? ^ 1.2). cf"^ shown by chain curves 
are overestimated compared to cf ^ and c^^. Dashed curves denoting c^^ [49] are 
plotted only for 0.8 < q < 1.2, because the FA is considered to be valid for a small 
|g — 1| [23j. The q dependence of c^^ is qualitatively different from those of the EA, 
lA and SA: c^^ is symmetric with respect to g = 1.0 whereas those in other three 
methods are monotonously increased with increasing q. 

The energy of electron systems at low temperatures may be calculated with the 
use of the generalized Sommerfeld expansion. By using Eqs. (11081) and (IllOp for 
Eq. (!73l) with 0(e) = ep(e), we obtain the energy given by 

E,iT) ^ E,iO) + C2,,{kBTfp{fi) + -, (118) 

from which the low-temperature electronic specific heat is given by 

C,{T) ^ 7,r + -, (119) 

with 

^ = ^, (120) 

7l C2,l 

71 = y4p(/"), (121) 

where 71 is the linear-T expansion coefficient for q = 1.0. 

The inset of Fig. [1] shows that the calculated energy Eg at low temperatures in 
the electron model is larger for a larger q, which is consistent with larger 7^ and C2,q 
for a larger q as shown in Fig. [TlT a). 

4.3 Low-temperature phonon specific heat 

We consider the phonon specific heat at low temperatures. By using Eqs. (!60l) and 
( 1651) . we obtain 



Ca ~ a,T^ + -, (122) 



with 



1 + 6(^-1) + -, (123) 



5 

where ai is the relevant coefficient for q = 1.0. 



127r''' 

«i = ( ^ ) NakB, (124) 
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The coefficients of low-temperature phonon specific heat aq in the lA, SA and 
FA are given by (for details, see Appendix B.3) 



ai (2-g)(3-2g)(4-3g)^ 
«f _ 1 



ai (2-g)(3-2g)(4-3g)(5-4g)' 



:i25) 



(126) 



= l + 0((g-lf), (127) 

where the 0{q — 1) contribution to a^^ is vanishing [19]. Equation (11251) shows 
that a^"^ agrees with a^^ within 0{q — 1) and that the a^"^ is related with c^^ as 
a^^/ai = cy/c4,i. 

Coefficients of Og/ai calculated by the four methods are plotted as a function 
of q in FigJTlTb). Squares denote the result of numerical calculation by the EA for 
T/Td = 0.01 (Fig. H]). The solid curve express a^^ which is in good agreement with 
the result of the EA for 1.0 < g < 1.2 but deviates from it at g > 1.2. Dashed and 
chain curves show aq calculated by the FA and SA, respectively. It is interesting 
that the result of the SA nearly coincides with that of the FA for 1.0 < g < 1.2, 
where both the results of the SA and FA are overestimated compared to the EA. 
The inset of Fig. H] shows that the energy Eq at low temperatures in the Debye 
model is larger for larger q, which is consistent with the g-dependence of aq shown 
in Fig. Mh). 



5 CONCLUDING REMARKS 

It is well known that in nonextensive classical statistics, the nonextensivity arises 
from the long-range interaction, long-time memory and a multifractal-like space-time 
[2]. The metastable state or quasi-stationary state is characterized by long-range in- 
teraction and/or fluctuations of intensive quantities [e.g., the inverse temperature) 
|10j . For example, in the long-range-interacting gravitating systems, the physical 
quantities are not extensive: the velocity distribution obeys the power law and the 
stable equilibrium state is lacking, which lead to negative specific heat p3] . The sit- 
uation is the same also in nonextensive quantum statistics. It has been reported that 
the observed black-body radiation may be explained by the nonextensivity of the 
order of |g — 1| ~ 10"'' — 10~^ which is attributed to the long-range Coulomb inter- 
action |21j. Memory effect and long-range interaction cannot be neglected in weakly 
non-ideal plasma of stellar core [M]- In addition to the large systems where the 
interactions may be truly long range, one should consider small systems where the 
range of the interactions is of the order of the system size. Small-size systems would 
not be extensive, and many similarities with the long-range case will be realized. 
Indeed, the negative specific heat is observed in 147 sodium clusters |65]. Mag- 
netic properties in nano- magnets may be different from those in large-size ones [SS] • 
Small drops of quantum fluids may undergo a Bose-Einstein condensation. Thanks 
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to recent development in the evaporation cooling technique, it becomes possible to 
study Bose-Einstein condensation in an extremely diluted fluid where the long-range 
interactions play essential roles in the condensate stability. Artificial sonic or optical 
black hole [SZl EB] represents an intrigue quantum catastrophic phenomenon. Only 
little is known about the thermodynamics of these quantum systems. Experimental 
and theoretical studies on these subjects deepen our understanding of basic quantum 
phenomena. 

To summarize, we have discussed the generalized distributions of g-BED and 
g-FDD in nonextensive quantum statistics based on the EA [50l [51] and lA. Results 
obtained are summarized as follows: 

(i) with increasing q above q = 1.0, the g-BED and g-FDD have long tails, while 
they have compact distributions with decreasing q from unity, 

(ii) the coefficients in the generalized Sommerfeld expansion, the linear-T coefficient 
of electronic specific heat and the coefficient of phonon specific heat are increased 
with increasing q above unity, whereas they are decreased with decreasing q below 
unity, 

(iii) the 0{q — 1) contributions in the EA agree with those in the AA based on the 
OLM-MEM [5] as well as the un-normalized MEM [3], and 

(iv) the generalized distributions given by simple expressions in the lA proposed in 
this study yield results in agreement with those obtained by the EA within 0{q — l) 
and high- and low-temperature limits. 

As for the item (iv), the g-BED and g-FDD in the lA are expected to be useful and 
to play important roles in the nonextensive quantum statistics. 
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A The {q - 1) EXPANSION IN THE UNNOR- 
MALIZED MEM 

Tsallis, Sa Barreto and Loh [21] developed the AA to investigate the nonextensivity 
in the observed black-body radiation, by using the un-normalized MEM |3j. We will 
show that the EA with the un-normalized MEM yields the result in agreement with 
the AA within 0{q — 1). Calculations of the g-BED and g-FDD for g ~ 1.0 are 
presented. 
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A.l Un- normalized MEM 



An application of the un-normalized MEM to the hamiltonian H yields the optimized 
density matrix given by [3] 



1 



Pi 



[l-[l-q)(3H]T-~ 



= Tr{[l-{l-q)m~^}- 
The expectation value of the operator O is given by 
0,{P) ^ (d), = Tr{p^O}, 

= i^Tr{[l-(l-g)/5iy]T^6}. 



(Al) 
(A2) 

(A3) 
(A4) 



A. 2 Exact approach 

With the use of the exact representations given by Eqs. ([5]) and (127|) . Eqs. (]A2I) 
and (1A4[) are expressed by 



oo / ^ 

G u: 



c 



0„ 





i 

2^ 
1 



1 



if t 



1 1 



G ( U- r + 



q JC 



q-l 
1 



-1, 



Zi{u) du for g > 1, (A5) 
Zi{-t) dt for g < 1, (A6) 
Zi{u)Oi{u) du for g > 1, (A7) 
Zi{-t)Oi{-t) dt for g < 1, (A8) 



1 



;i - q)p 



with 



Oi{u) 



Tr{e-''^ 6} 



(A9) 
(AlO) 



where C denotes the Hankel contour, and G{u; a, b) and H{t; a, b) are given by Eqs. 
(fni) and (I5T]) . respectively. In order to evaluate Eqs. flASp - flASp . we expand their 
integrands around u = P and — t = /5 as is made in Sec. 2.3. By using Eqs. (!22l) . 
(123]), (BSD and (1361), we obtain 



d^Z, 



+ ■■ 



On 



^q 



Oi + (g - I)f3-^{Z,0,) + i(g - 1)(3'-^{Z^0^) 



(All) 
• (A12) 
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By using the relations given by 
dZi 



d(3 
dOi 

'W 



1) 



(/J)i(0)i-(i70)i, 



we finally obtain the 0{q — 1) expansion of Og given by 

~ 0, + {l~q)(^Oi\nZ, + P{Hd), + ^P^[{H^)iOi-{H^d)i]^+--, 

(A13) 

which agrees with Eq. (7) of Ref. [21] derived by the AA. 
(1) g-BED 

In order to calculate the g-BED, we consider O = fik with the hamiltonian for 
bosons given by 

^ = ^(^k - /i) nfc, (A14) 

k 

where and stand for the number operator and the energy of the state k. We 
obtain 

{nk)i = -^ = h{ek) = fu [x = (3{ek-fi)] (A15) 

(rifcif)i = {ek-fi)e^f! + fiE,, (A16) 
= ^2 + ^2 + ^3, (A17) 
{fikH'), = 2{tk-^ffl-2{tu-^)flE^ + h{El + E2 + E^), (A18) 

with 

El = ^(efc - 

k 

E2 = 5^(efe-/i)Vi, 

k 
k 

Substituting Eqs. flXTSD - flXTSj) to Eq. fIXTSD . we obtain 

/, ^ /i + (l-g)(/ilnZi+/3[(e,-/i)eVi' + /ii5i]) 

- ^i^^^[(6,-/i)V(e^' + l)A3 + 2(6,-^)eVi'i?i]+--. (A19) 

24 



Tsallis et. al. [21] employed a one- component boson hamiltonian given by 



H 



which yields 



{nk)i 
{nH)i 



1 



tujj n = eh. 



{x = (3e) 



A substitution of Eqs. (IMT|) - (IA24D to Eq. (IM3|) leads to 



/i + (1 - q) 



h InZi + x{e^ + l)fl - -x2e^(e^ + 3)/f 



+ 



(A20) 



(A21) 

(A22) 
(A23) 
(A24) 



(A25) 



which is different from Eq. (]A19P with ^ = because of the difference in the adopted 
hamiltonians given by Eqs. (]A14I) and (]A2ip . 



(2) g-FDD 

We consider O = hk with the hamiltonian for fermions given by 

H = ^(e/c - ft) hk, 



(A26) 



which leads to 

{hkH)i 
{hkH^)i 



= A(efc) = A, [a; = - /i)] (A27) 

(efc-/i)/i(l-/i) + /iEi, (A28) 

El + E2-E3, (A29) 

(efc - - /i)(l - 2/1) + 2(efc - - h)E, 



+ fi{El + E2-Es). (A30) 
Substituting Eqs. (lA27j) -( lA30ll to Eq. (lATl) . we obtain 

/, ^ /, + (i_g)(/^lnZi + /?[(efc-/i)/i(l-/i) + /iEi]) 

- ^i^^ [(6, - - - 2/1) + 2(6, - - + ■ • . 

(A31) 



When assuming a one-component fermion hamiltonian given by 

H = {ek- fi)hk, 



(A32) 
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we obtain 



{nkH)i 



1 



= f- 

e^' + 1 
(efc - 
(efc 



1) 



Substituting Eqs. flA33l) - flA36|) to Eq. flATSl) . we obtain 

1 



/i + (1 - g) 



(A33) 

(A34) 
(A35) 
(A36) 



.(A37) 



The difference between Eqs. (IA3ip and flA37l) is due to the difference in the adopted 
hamiltonians given by Eqs. flA26l) and flA32|) . It is noted that the (g - 1) expansion 
of g-FDD in the FA is given by 



FA 



fl 



/52(e-/i)V(^-^)/2 + 



(A38) 



whose 0{q — 1) term corresponds to the last term of Eq. (]A37p derived by the 
un-normahzed MEM. This is due to the fact that to adopt the one-component 
hamiltonian given by Eq. (lA32p means to use the factorization approximation from 
the beginning. 

Equation flAlQP for g-BED and Eq. flA3ip for g-FDD are expressed in a unified 
way as 



fg ^ fl 

+ (1-9) 



+ -, 
(A39) 



where /i = l/(e^=Fl). We note that the 0{q — l) term of the generahzed distribution 
in Eq. fl65|) derived by the OLM-MEM corresponds to the last term in the bracket 
ofEq. flA39l) . 



B SUPPLEMENT TO THE INTERPOLATION 
APPROXIMATION 

B.l Analytic expressions of g-FDD for \/3{e — fi)\ <^ 1 

We may obtain an expression of the g-FDD for small x [= (3{e — yu)] with the use of 
an expansion for /i(e/3) given by 

= 2+5Z<i^" for|x|<l. (Bl) 

n=l 
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Substituting Eq. ([BT]) to Eqs. ([76]) and ([77]), and employing Eq. ([5]) and ([27]), we 
obtain 

^ oo 

Cie,P) = i^ + E^-?^" for|x|<l, (B2) 



n=l 



with 



(g-l)"r(^ + 1 + n) 

dn,q = dn,i ^ for 1 < g < 3, (B3) 

(1 -g)"r(^ + 1) 

= <i p. , ^ , forO<g<l, (B4) 

= q dn,i for n = 1, (B5) 

= q{2q - 1) for n = 2, (B6) 

= g(2g-l)(3g-2Xi forn = 3, (B7) 

where dn,i = (1/n!) (9"/i(e, /5)/(9x" at x = 0: cii,i = -1/4, c/2,1 = 0, 4,1 = 1/48, 
etc.. Equations ([B2])-([Mj) lead to 



3) . + forW<l. (B8) 



B.2 Generalized Sommerfeld expansion in the I A 

In the case of g > 1.0, Eq. (|6T1) yields 

Substituting Eq. ( ]B9p to Eq. ( 1109P and changing the order of integrations for e and 
u, we obtain 

At low temperatures, Eq. (IBlOp reduces to 

2(l-2i-")CH q \_ 

r(^ + l-r2) 

= (g-Vr(^ + i) ^""''^ 

= for odd n. (B13) 



The ratio of Cn,q/cn,i is given by 



c r(^ + l-n) 



(g - i)«r(^ + 1) 
1 

2^ 

1 



(2-g)(3-2g)(4-3g) 
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for even n, (B14) 

for n = 2, (B15) 
for n = 4. (B16) 



In the case of g < 1.0, Eqs. ([62]) and (fT09|) yield 

2(l -2^-")CH I [ ^(^^ 
[l-qY 27rJc \'l-q 



cn,, = jr^ ^ I H(t,-r^A] i-t)- dt, (B17) 



leading to 



r(T^ + i) 

for odd n, (B19) 



= , for even n, (B20) 

Cn,i (l-g)«r(^ + l + n) ^ ' 

1 



2-q 

1 



for n = 2, (B21) 

r forn = 4. (B22) 

(2-g)(3-2g)(4-3g) ^ ' 

Equation flB20D for g < 1.0 is the same as Eq. (]B14p for g > 1.0 if we employ the 
reflection formula of the gamma function: 

rrz)rri-^) = 



Sm TTZ 



B.3 The low-temperature phonon specific heat in the I A 

In the case of g > 1.0, Eqs. fl60l) and fl76l) yield 



oo 







Ca ^ ksP^ / G{u; — ^— , 1 / r rJ.U. du du 



g - 1 ' J Jo [e(g-i)/3S^« _ 1] 



" ^ / G i u\ — ^— , 1 I M ^ du [ -, —Trdx, 



(B23) 



rj-] \ 3 



«.(^) , (B24) 



with 



r(^-3) 

aq = ai for 1 < g < 3 , (B25) 

where To (= huJo/kB) stands for the Debye temperature and ai is the coefficient 
of the low-temperature specific heat for g = 1.0. 

In the case of g < 1.0, a similar analysis with the use of Eqs. fl60l) and (I77|) leads 

to 



(B26) 



[l-qY\TDj \2TiJJc V y io (e^--l)' 
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from which we obtain 



an 



r(A + i) 



(l-g)4r(^ + 5) 
Equations (IBTHI) . (!B22|) . (!B251) and (!B28l) yield 



for < g < 1. 



(B28) 



C4,, 



(2-g)(3-2g)(4-3g) c,,_ 



for < g < 4/3. (B29) 
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Table 1: Generalized distributions in the limits of g — > 1, T — and (3^0 



method 




T ^ (FDD) 


p ^ 


EA" 


fi + {q-l) 




+ ■■ 


e(/i-e) 




lA^ 


fi + {q-l) 




+ ■■ 


e(/i-e) 




FA'^ 


/i-i(g-l)/3(e-/.)2^ + - 


0(/i-e) 




SA"^ 




0(/i-e) 


^-/3(f-M) 



/i = l/(e'^'''^~'^'' =F 1): 0(x), the Heaviside function: e^, g-exponential function. 

" the exact approach (the present study) 

^ the interpolation approximation (the present study) 

the factorization approximation ^24] 

the superstatiscal approximation [39] 
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Table 2: 0{q — 1) contributions to c„^g (ra = 1 — 4) of the generalized Sommerfeld 
expansion coefficients 



method 




C2,q 


C3,q 


C4,g 












+ 6(9-1)] 


lA^ 





^[1 + (?-!)] 





S[l + 6(9-1)] 


FA" 




f[l + 0((g-in] 




g[l + 0((g-in] 


SA*^ 





T[i + 3(g-l)] 





||;[l + 10(g-l)] 



" the exact approach (the present study) 
the interpolation approximation (the present study) 
the factorization approximation 
the superstatiscal approximation 
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Figure 1: (Color online) The temperature dependence of Eg of the electron model 
for q — 1.0 (dashed curves), q — 1.1 (chain curves), q — 1.2 (dotted curves) and 
q — 1.3 (sohd curves), the inset showing the enlarged plot for ksT/W < 0.1. 



Figure 2: (Color online) The e dependence of the g-FDD of fq{e) for q = 0.8 (solid 
curves), q — 0.9 (dotted curves), q = 1.0 (dashed curves), q = 1.2 (double-chain 
curves), q — 1.5 (bold solid curves) and g = 1.8 (chain curves) with (a) the linear 
and (b) logarithmic ordinates, the results for q > 1.0 and q < 1.0 being calculated 
by the EA and lA, respectively {kBT/W = 0.1). 



Figure 3: (Color onhne) The e dependence of the derivative of g-FDD, —dfq{e)/de, 
for q — 0.8 (the sohd curve), q — 0.9 (the dotted curve), q — 1.0 (the dashed curve), 
q — 1.2 (the double-chain curve), q — 1.5 (the bold solid curve) and q = 1.8 (the 
chain curve) with the logarithmic ordinate, the results for q > 1.0 and q < 1.0 being 
calculated by the EA and lA, respectively {ksT/W = 0.1). 



Figure 4: (Color online) The temperature dependence of Eg of the Debye phonon 
model for g = 1.0 (dashed curves), q = 1.1 (chain curves), q = 1.2 (dotted curves) 
and q = 1.3 (solid curves), the inset showing the enlarged plot for T/Tp < 0.5. 



Figure 5: (Color online) The e dependence of the g-BED of fq{e) for q — 0.8 (the 
solid curve), q — 0.9 (the dotted curve), q — 1.0 (the dashed curve), g = 1.1 (the 
chain curve), g = 1.2 (the double-chain curve), g = 1.5 (the bold solid curves) and 
g = 1.8 (the thin solid curve) with the logarithmic ordinate, the results for g > 1.0 
and g < 1.0 being calculated by the EA and lA, respectively (T/Td = 0.01). 



Figure 6: (Color online) The e dependence of the g-FDD of fq{e) calculated by the 
EA for g = 1.0 (dashed curves), g = 1.2 (chain curves), g = 1.5 (dotted curves) and 
g = 1.8 (solid curves) with the logarithmic ordinate, the inset showing the ratio of 
A = /<;'"(e)//<f-^(e) {kBT/W = 1.0). 
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Figure 7: (Color online) The e dependence of the g-BED of /g(e) calculated by 
the EA for g = 1.0 (dashed curves), g = 1.1 (double-chain curves) q = 1.2 (chain 
curves), q = 1.5 (dotted curves) and q = 1.8 (solid curves) with the logarithmic 
ordinate, the inset showing the ratio of A = fq^{^)/ fq^{^) {ksT /W = 0.1). 



Figure 8: (Color online) The e dependence of the g-BED of fq{e) for g = 1.1 and 
1.2 calculated by the EA (solid curves), FA (chain curves) and SA (dotted curves) 
with the logarithmic ordinate, /i(e) for g = 1.0 being plotted by the dashed curve 
for a comparison (T/Td = 0.01). 



Figure 9: (Color online) The e dependence of the g-FDD of fq{e) for g = 1.1 and 
1.2 calculated by the EA (the solid curve), FA (the chain curve) and SA (the dotted 
curve) with the logarithmic ordinate, /i(e) for g = 1.0 being plotted by the dashed 
curve for a comparison [ksT /W = 0.1). 



Figure 10: (Color online) The e dependences of (a) the g-FDDs of fq{e) and (b) 
its derivative of —dfq{e)/de calculated by the lA for g = 0.9 (sohd curves) and 1.1 
(bold solid curves), and those calculated by the FA for g = 0.9 (dashed curves) and 
1.1 (bold dashed curves), results for g = 1.0 being plotted by chain curves for a 
comparison. 



Figure 11: (Color onhne) (a) The g dependence of c„,q/c„_i for n = 2 and 4 of 
the generalized Sommerfeld expansion coefficients [Eq. fll08p ] with the g-FDD, and 
(b) the g dependence of Og/ai of the coefficients in the low-temperature phonon 
specific heat with the g-BED, calculated by the EA (circles and squares), lA (solid 
curves), FA (dashed curves) |19] and SA (chain curves): the result of the SA is 
indistinguishable from that of the FA in (b) (see text). 
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